Stimulated scattering and lasing of intersubband cavity polaritons 
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We present a microscopic theory describing the stimulated scattering of intersubband polariton 
excitations in a microcavity-embedded two-dimensional electron gas. In particular, we consider the 
polariton scattering induced by the spontaneous emission of optical phonons. Our theory demon- 
strates the possibility of final-state stimulation for the scattering of such composite excitations, 
accounting for the deviations from perfect bosonicity occurring at high excitation densities. By 
using realistic parameters for a GaAs semiconductor system, we predict how to achieve a quantum 
degenerate regime, leading to ultralow threshold lasing without population inversion. 



The scattering of bosons from an initial to a final state 
can be stimulated, i.e., enhanced by the occupation of 
the final state. This remarkable property is in stark 
contrast with the behavior of fcrmions, such as elec- 
trons, whose scattering is Pauli blocked by final state 
occupation. In low-energy matter there are no elemen- 
tary bosons, yet composite particles acting like bosons 
can be obtained when an even number of fermions are 
bound together, such as atoms containing an even total 
number of nucleons plus electrons. In condensed mat- 
ter systems [13], the attractive interaction between two 
electrons can give rise to bosonic particles. Examples 
are Cooper pairs of electrons in metallic superconductors 
or Coulomb bound electron-hole pairs (excitons) in semi- 
conductors. The strong coupling of an exciton with a mi- 
crocavity photon produces the so-called exciton-polariton 
states, whose very small mass (inherited from the photon 
component) favors quantum degeneracy and the onset of 
stimulated scattering, responsible for exciton-polariton 
lasers [2] emitting in the near infrared. 

Recently a novel kind of cavity polariton excitations 
has been discovered in a microcavity-embedded two- 
dimensional electron gas, [|| and an intense research ac- 
tivity is currently expanding 0, H, S 0, 11 S E, El E3 • 
These light-matter excitations are the result of the strong 
coupling between a microcavity photon mode and the 
transition between two conduction subbands of the doped 
quantum well system. A intersubband excitation has 
a well definite resonance frequency simply because the 
quantum well conduction subbands have parallel energy- 
momentum dispersions (see Fig. [T]). In contrast to 
Cooper pairs or excitons, intersubband excitations do 
not correspond to any bound state originated from an 
attractive fermion-fermion interaction. This explains the 
remarkable robustness of intersubband cavity polaritons 
even at room temperature and the possibility of accu- 
rately tailoring their properties just by tuning the size of 
the quantum well or the density of the two-dimensional 



electron gas in the fundamental subband. Even if inter- 
subband excitations do not correspond to any bound elec- 
tronic states, a priori this fact does not preclude the oc- 
currence of stimulated scattering. Yet, to the best of our 
knowledge, a theory for the stimulation phenomena in- 
volving these composite excitations is still missing. Apart 
from the fundamental interest, a comprehensive under- 
standing of the stimulated scattering of intersubband 
cavity polaritons could pave the way to the remarkable 
realization of ultraefficient semiconductor lasers emitting 
in the mid and far infrared. 

In this Letter, we present a microscopic theory of the 
stimulated scattering of intersubband cavity polariton 
excitations of a two-dimensional electron gas. In par- 
ticular, we will consider the polariton scattering induced 
by the coupling with optical phonons, which is typically 
the most important interaction affecting semiconductor 
intersubband transitions, while Coulomb interactions are 
known to produce only moderate renormalization effects 
[T3 |. Starting from the fermionic Hamiltonian for the 
quantum well electronic system and by using an exact 
iterative commutation procedure, we are able to deter- 
mine the phonon-induced polariton scattering for an ar- 
bitrary number of excitations in the initial and final in- 
tersubband cavity polariton modes. Our results indeed 
prove the possibility of final-state stimulation of the in- 
tersubband cavity polariton scattering. Our theory also 
provides exactly the deviations from perfect bosonicity, 
occurring at high excitation densities. We apply our re- 
sults to the case of a GaAs system with realistic losses 
and consider the case of intersubband cavity polariton 
lasing under resonant optical pumping. 

We consider the Hamiltonian H = Hi m + H p hon 
where Hi m is the light-matter term for the microcavity- 
embedded quantum well system, while H p h on describes 
the dynamics of bulk optical phonons and their the cou- 
pling to the quantum well electrons. Namely: 
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kj = l,2 q k,q 
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where ct k , a q and <i q ?a are the creation operators re- 
spectively for an electron in the quantum well conduction 
subband j with in-plane wave vector k, a cavity photon 
with in-plane wave vector q and an optical phonon with 
three-dimensional wave vector (q, q z ) . Their respective 
energies are tiVj k, ^>c,q and f)iOLo,q,q z (the wave vector 
dependence of the optical phonon energy is negligible) 
and their phases are chosen in order to make the cou- 
pling coefficients x(q) an( i CV?',q,<? z real. Being all the 
interactions spin conserving, we omit the spin degree of 
freedom for the electrons. The photon polarization is 
meant to be Transverse Magnetic (TM) in accordance 
with the selection rules of quantum well intersubband 
transitions. Note that, neglecting the conduction band 
non-parabolicity, the second subband dispersion is such 
that ujf k = Lof k + UJ12, as depicted in the inset of Fig.l. 
Moreover, for typical photonic wave vectors q, we can 
safely approximate w^ k+q | ~ ^fk- The light-matter 
Hamiltonian Hi m is then diagonalized by introducing the 
polariton creation operators 

Pr),q = a tJ,q a q Prj,qbg (2) 

where r\ = {LP, UP} denotes the polariton branch index, 
a^q and /S^q are real Hopficld coefficients describing the 
light and matter component respectively, while hhj^.q are 
their corresponding energies (see Fig. 1). 6 q is given by 
the expression 

h \ = 7^H C 2,k+q C l,k (3) 

with N the total number of electrons in the doped quan- 
tum well. In the ground state, the two-dimensional elec- 
tron gas fills the fundamental subband for k < kp, being 
kp the Fermi wave vector. & q creates a bright intersub- 
band excitation with in-plane wave vector q, obtained 
when the two-dimensional electron gas absorbs one cav- 
ity photon, as it can be deduced from the light-matter 
coupling in Eq. |T]). 

We are interested in calculating the polariton scatter- 
ing rate induced by the emission of an optical phonon 
from an initial polariton 'pump' mode (branch rj' and in- 
plane wave vector q') to a final 'signal' mode (branch r\ 
and in-plane wave vector q) . This kind of process is pic- 
tured in Fig. [I] for the case rf = UP and r\ = LP. It is 
clear that in order to have a sizeable polariton-phonon in- 
teraction, both the initial and final polariton modes must 



H 0.8 



V Pump 

In-plane wave vector At 








\ \ y/ Signal 




/ LP 





1.6 1 " 1 1 1 

0.5 1 1.5 2 



In-plane wave vector q (in units of q res ) 

FIG. 1: A typical energy dispersion (in units of the inter- 
subband transition energy hwi2) of intersubband cavity po- 
laritons versus in-plane wave vector (in units of the resonant 
wave vector g res ). Due to the interaction with bulk optical 
phonons, a polariton initially pumped in the upper polari- 
ton (UP) branch can scatter into a final state (signal mode) 
in the lower polariton (LP) branch by emitting an optical 
phonon with energy huiLO (36 meV for GaAs). The consid- 
ered modes have Hopheld coefficients Pup,^' = Plp,^ = 0.5. 
The dashed lines indicate the same kind of scattering process 
by changing the in-plane momentum of the initial state along 
the upper polariton branch. Inset: the energy dispersion of 
the quantum well electronic conduction subbands versus elec- 
tron wave vector k. In the ground state, a dense electron gas 
populates the first conduction subband. 



have significant electronic components, quantified respec- 
tively by \Pup,q'\ 2 and \0LP,q\ 2 - At the same time, in 
order to have a good coupling to the extracavity electro- 
magnetic field (required for pumping and detection) also 
the photonic components |o;f/p,q'| 2 an d I^Lpql 2 need to 
be significant. These conditions can be simply met when 
the polariton energy splitting 2hx(q res )^/~N at the reso- 
nant wavector q res (such as uj c (q res ) = LO12) is a non neg- 
ligeable fraction of the optical phonon energy (36meV for 
GaAs). This situation is already realized in recent micro- 
cavity samples [a, |9|, with mid-infrared intersubband 
transition frequencies. 

If we wish to investigate the occurrence of stimulated 
scattering, we need to evaluate the scattering rates for 
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arbitrary occupation numbers m and n of respectively 
the initial and final polariton modes. The emission of 
an optical phonon can induce the scattering of one po- 
lariton from the pump to the signal mode, leading to 
a transition from the state pJ^q'Pt^F) to the state 

^q'-q 9z^'"q~ lj 4,q +1 l^ where \F) is the ground state of 
the system (the electronic ground state times the photon 
and phonon vacuum). Therefore, we need to consider the 
squared normalized matrix element ft 2 |V^| 2 given by 



K^KWCq'^' 



-q,<?. 



»7,qi 



•(4) 



In order to evaluate Eq. 0$, we have to exploit the ex- 
pression of the polariton operators in Eq. ([2]) in terms of 
the cavity photon and intersubband excitation operators. 
To evaluate the matrix elements, we need to commute the 
destruction operators multiple times to the right side and 
exploit the annihilation identity a^\F >= bq\F >= 0. 
The cavity photons are elementary bosons obeying the 
standard commutation rule [a q ,a q/ ] = <5 q , q '- Instead, 
intersubband excitation operators arc not elementary 
bosons and satisfy modified commutation rules. We have 
found that 



b a> b U = <W _ D q.q'' 



q> q 
D, 



(5) 



f 



q,q 



2.^ C l,k C l,k+q-q' - c 2,k+q' c 2,k+q, 
|k|<k F 



where -D q , q ' is the operator describing the deviation from 
the behavior of elementary bosons, originally introduced 
in the context of excitonic composite bosons [15j. By iter- 
ation, we have found the following commutation relations 



[-t-Aj.q' , U q » J — N U q"+q'-q U q" > 



(6) 



[°q'°q' J '' t0 q' l°q,q' u 1,q.' ) N 2q'-q U q' 

[b^,b^}=m{5^ > -D^)h l 



fm-2 



m — 1 



m(rn— 1) j 
TV 
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Due to the fact that typical photonic wave vectors q are 
much smaller (at least two orders of magnitude) than the 
Fermi wave vector ftp, we have D^q>\F) ~ with cor- 
rections of the order of |q — q'lfcp- due to the electrons 
occupying the edge of the Fermi sphere in the ground 
state. Neglecting these corrections is the only approx- 
imation we will assume, which becomes exact in the 
limit |q — q'|fci? — » 0. Exploiting Eq. (O some algebra 
shows that the unnormalized polaritonic matrix element 
(F\^P^d^-^H ph ^y^pll\F) is given by 



(» + l)!m!A, l «A' J ,'(C2 2 , q - q ' ig ,-Cr lliq -^, 9 .) £ (") ("\ ^ l^l aI IA J ,q| a(n "°K',q'l a ^/V,^ 2 ^ -1 "^'/^) 



J=0,...,n 
/i=0,...,m-l 

where /" = + /C™_ ) ~ L 1 ^~ 1 and the quantity /C^ s r is defined by the relation 

n\m\K^ r = (F\b%b™,bl%rt Q \F) (8) 
with Q = q(n — s) + q'(m — r). Analogously for the normalization factors in Eq. ((¥]), we find 

W,^* = E (^(^K^l/W^ -0 !^ (9) 

/i— 0,...,m 



It is clear that JC r ^ r oc (5„+ m , r +s+i and that = 
m-i an d r = ^rTfT- Using the commutation rela- 
tions in Eq. ^ we obtain the the recurrence relation 

v~n s £ x yTi—l,n—l . c c i^n.n — l 

^m,r ~ Orn.rOn.s+ll^m^m-i + "m.,r+l"rc,s"''77l— l,m— 1 

that allows us to numerically evaluate the IC™'^ for arbi- 
trary and realistic occupation numbers. Hence, Eq. Q 



can be rewritten in the simpler form 

|C| 2 = (n + l)m^|A J> «A',q'(C3a,| q -q'|, fl . - Cu,| q -q'|, to )| 2 

where B 1 ^ is a bosonicity factor depending on the coef- 
ficients /Cm S r . Its expression is cumbersome, but it can 
be simply obtained putting together Eqs. (j4|), and 
dH) . Such a quantity depends on the Hopfield coefficients 
and on excitation numbers m and n normalized to the 
total number of electrons N in the ground state. In the 
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neglecting the LO-phonon dispersion, we obtain 
T sc = K n + 1 ) B m\Pr,,q\ l/V,q'l — -F 7i 1 11 ' 
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FIG. 2: Spontaneous scattering rate r™ ,n=0 for the process 
depicted in Fig. 1 versus the pump polariton density m/S. 
The electron density in the ground state is N/S = 10 12 cm -2 . 
In the considered range of excitation densities, m/N < 0.25, 
i.e., much smaller than the onset of electronic population in- 
version. Other GaAs parameters are given in the text. Inset: 
the solid line represents the bosonicity factor B™ =0 versus 
m/N for the pump and signal polariton modes considered in 
Fig. 1. For elementary bosons BJ, is always 1. The dashed 
line is the same quantity for pure matter excitations. For 
m/N -C 1, deviations from perfect bosonicity are negligible. 
For n <^ m (signal much smaller than pump), B 7 ^ ~ 
and the stimulated scattering rate is T™'™ ~ (1 + n)F™' n=0 . 



inset of Fig. 2, we report B^ versus m/N obtained by a 
numerical evaluation of our exact recursive relations. For 
normalized excitation densities m ^ n smaller than 0.1, we 
find that £?" is well approximated by the formula 



B n m * 1 - C- 



N 



(10) 



where £ depends on the Hopfield coefficients of the po- 
lariton modes and varies from for pure photonic exci- 
tations to 1 for pure matter ones. Hence, we see that 
final-state simulation, due to the enhancement factor 
(n + 1), can be obtained when the density of polaritons 
(m + n)/S is much smaller than the density N/S of the 
two-dimensional electron gas, not a stringent condition 
as shown later. Indeed, the bosonicity of intersubband 
cavity polaritons can be controlled by the changing the 
density of the two-dimensional electron gas, which can 
be tuned by an electric gate [|| . Using Fermi golden rule 
we find the following formula for the scattering rate 

r™'"=2^ [ dw|CrA,-q',9.,^Kq-^',q'+w), 

where A q _ q /^. w is the spectral function of the optical 
phonon. Using a Lorentzian shape of width Ylq and 
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where S is the sample surface, Lqw the quantum 
well length and F a a form factor (depending on a = 
Lqw\i — q'l) describing the overlap between the conduc- 
tion subband envelope functions and the optical phonon 
wavefunction [ljj] . For typical quantum well widths and 
photonic wave vectors, a <C 1. In the case of a quan- 
tum well with infinite barriers, F a ~Q ~ 0.1. For GaAs 
optical phonons, the ratio « 100 [17]. In Fig. H we 
report the calculation of the spontaneous scattering rate 
r™'° (i.e., n — 0, unoccupied final state) for the process 
showed in Fig. Q]for a GaAs system with hu>\2 = 150meV 
(mid infrared) and N/S — 10 12 cm~ 2 . The scattering 
rate into the final signal mode for pump excitation den- 
sities around 10 11 cm -2 is comparable to the typical loss 
rate of intersubband cavity polariton modes. Hence, for 
these pump excitation densities one expects to enter the 
regime of stimulated scattering. 

Neglecting the pump depletion (relevant only above an 
eventual stimulation threshold), we can write two rate 
equations for the signal and pump populations, namely 



dn 
~dt 



pm,n p 

— sc loss 



dm 
~~dt 



loss^ 



UP,q' 



(12) 



where T 



and r!„„„ are the loss rates of the 



sig- 



nal and pump modes, A the polariton absorption co- 
efficient at the pump frequency and I pU mp the optical 
pump intensity. From the steady-state solution for n, 
we can calculate the threshold pump density mthr/S 
to have a lasing instability. For n -C m, ~ B® n 
and r™<" ~ (1 + n)T™>°. The threshold pump po- 
lariton density mthr/S is then given by the equation 
T™ thr '° — Tioss ■ The steady-state solution for m gives the 
threshold pumping intensity versus the polariton thresh- 
old density, namely I^ mp = Vl -^ h ^ up ^ m thr /S. For a 
realistic Ti oss — 5 ps _1 , we obtain a threshold density 
for the pump mode of 1.1 x 10 n cm~ 2 , i.e. m/N = 0.11, 
as indicated in Fig. 2. With a polariton absorption 



coefficient A = 0.4 11], this gives an impressively low 
threshold pump intensity of 3.5 x 10 4 IU/cm 2 . This is 
approximately 2 orders of magnitude smaller of what re- 
quired to achieve electron population inversion in the two 
subbandsfljj]. 

In conclusion, we have derived a theory for the stim- 
ulated scattering of intersubband cavity polariton exci- 
tations of a dense two-dimensional electron gas. The 
intersubband cavity polariton excitations are composite 
bosons arising from the strong light-matter coupling and 
are not associated to any bound electronic states. We 
have shown exactly how the bosonicity of these excita- 
tions is controlled by density of the two-dimensional elec- 
tron gas in the ground state. Our exact results show that 
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ultralow threshold intersubband cavity polariton lasing 
without population inversion can be achieved under op- 
tical pumping by exploiting the efficient interaction with 
optical phonons. The present theory could pave the 
way to the experimental demonstration of fundamental 
quantum degeneracy phenomena and unconventional las- 
ing devices based on composite bosons with controllable 
properties and interactions. 

We wish to thank S. Barbieri, R. Colombelli, F. Julicn 
and C. Sirtori for discussions. 



[1] W. Ketterle, Rev. Mod. Phys. 74, 1131 (2002). 

[2] Le Si Dang et al, Phys. Rev. Lett. 81, 3920 (1998); J. 

Kasprzak et al, Nature 443, 409 (2006). 
[3] D. Dini et al, Phys. Rev. Lett. 90, 116401 (2003). 
[4] C. Ciuti, G. Bastard, I. Carusotto, Phys. Rev. B 72, 

115303 (2005). 

[5] R. Colombelli et al, Semicond. Sci. Technol. 20, 985 
(2005). 

[6] A. A. Anappara et al, Appl. Phys. Lett. 87, 051105 



(2005); A. A. Anappara et al, Appl. Phys. Lett. 89, 
171109 (2006). 

[7] C. Ciuti, I. Carusotto, Phys. Rev. A 74, 033811 (2006). 
[8] L. Sapienza et al, Appl. Phys. Lett. 90, 201101 (2007). 
[9] A. A. Anappara et al, Appl. Phys. Lett. 91, 231118 
(2007). 

[10] S. De Liberato, C. Ciuti, I. Carusotto, Phys. Rev. Lett. 
98, 103602 (2007). 

[11] L. Sapienza et al, Phys. Rev. Lett. 100, 136806 (2008). 

[12] S. De Liberato, C. Ci uti, Phys. Rev. B 77, 155321 (2008) 
and larXiv:0802.409Tl 

[13] A. J. Leggett, Quantum Liquids: Bose Condensation and 
Cooper Pairing in Condensed-Matter Systems, (Oxford 
University Press, 2006). 

[14] D. E. Nikonov et al., Phys. Rev. Lett. 79, 4633 (1997). 

[15] M. Combescot, O. Betbeder-Matibet, R. Combescot, 
Phys. Rev. B 75, 174305 (2007); O. Betbeder-Matibet, 
M. Combescot, Eur. Phys. J. B 27, 505 (2002) and ref- 
erences therein. 

[16] R. Ferreira, G. Bastard, Phys. Rev. B 40, 5802 (1989). 

[17] P. Y. Yu, M. Cardona, Fundamentals of Semiconductors, 
(Springer, 2001). 

[18] O. Gautier-Lafaye et al, Appl. Phys. Lett. 71, 3619 
(1997). 



